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Abstract 

Higgs-boson pair production is well known being capable to probe the trilinear self-coupling of 
the Higgs boson, which is one of the important ingredients of the Higgs sector itself. Pair pro- 
duction then depends on the top-quark Yukawa coupling ’ , Higgs trilinear coupling X^h, and 
a possible dim-5 contact-type ttHH coupling , which may appear in some higher representa¬ 
tions of the Higgs sector. We take into account the possibility that the top-Yukawa and the ttHH 
couplings involved can be CP violating. We calculate the cross sections and the interference terms 
as coefficients of the square or the 4th power of each coupling [g^ ’ , Xsh, Qtt ) various stages 
of cuts, such that the desired cross section under various cuts can be obtained by simply inputing 
the couplings. We employ the HH ^ 7766 decay mode of the Higgs-boson pair to investigate the 
possibility of disentangle the triangle diagram from the box digram so as to have a clean probe of 
the trilinear coupling at the LHC. We found that the angular separation between the b and b and 
that between the two photons is useful. We obtain the sensitivity reach of each pair of couplings 
at the 14 TeV LHC and the future 100 TeV pp machine. Finally, we also comment on using the 
hhT^T~ decay mode in Appendix. 
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I. INTRODUCTION 


A boson was discovered at the Large Hadron Collider (LHC) [U [2] . After almost all the 
Run I data were analyzed, the measured properties of the new particle are best described 
by the standard-model (SM) Higgs boson 015) which was proposed in 1960s O- The most 
constrained is the gauge-Higgs coupling = gnww = 0-94 which is very close to the 
SM value [6]- On the other hand, the relevant top- and bottom-Yukawa couplings are not 
determined as precisely as by the current data. Nevertheless, they are within 30 — 40% 
of the SM values [6] . 

Until now there is no information at all about the self-couplings of the Higgs boson, which 
emerges from the inner dynamics of the Higgs sector. For example, the trilinear couplings 
from the SM, two-Higgs doublet models (2HDM), and MSSM are very different from one 
another. Thus, investigations of the trilinear coupling will shed lights on the dynamics of 
the Higgs sector. One of the best probes is Higgs-boson-pair production at the LHC. There 
have been a large number of works of Higgs-pair production in the SM [MO], in model- 
independent fashion [TTllTT] . and in special models beyond the SM [I5] and in SUSY [T6] . 
In the SM, Higgs-pair production receives contributions from two entangled sources, the 
triangle and box diagrams. The triangle diagram involves the Higgs self-trilinear coupling 
and the top-Yukawa coupling while the box diagram involves only the top-Yukawa coupling. 
In order to probe the effects of the Higgs trilinear coupling, we have to disentangle the 
triangle diagram from the box diagram. We anticipate that the triangle diagram, which 
contains an s-channel Higgs propagator, does not increase as much as the box diagram as 
the center-of-mass energy a/S = Mhh increases. Therefore, the box diagram tends to give 
more energetic Higgs-boson pairs than the triangle diagram. Thus, the opening angle in the 
decay products of each Higgs boson can be used to isolate the triangle-diagram contribution. 
Indeed, we found that the angular separation and AR^b between the decay products 

of the Higgs-boson pair are very useful variables to disentangle the two sources. 

Here we also entertain the possibility of a dimension-5 operator ttHH, which can arise 
from a number of extended Higgs models, including composite Higgs models or some general 
2HDM’s. For example, in a general 2HDM we can have a diagram with a {titRip) vertex 
and a {ipHH) vertex connected by the heavy ip. When the heavy ip is integrated out, we are 
left with the contact diagram titRlIH. The anomalous ttHH coupling can contribute to 
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Higgs-pair production via a triangle diagram. This triangle diagram is similar to the triangle 
diagram with the trilinear Higgs self coupling, except that it does not have the s-channel 
Higgs propagator. We shall show that the new contact diagram will give terms that can be 
combined with the terms of the triangle diagram, as in Eq. Q. We note that the kinematic 
behavior of the triangle diagram induced by the dim-5 ttHH contact interaction is different 
from that induced by the trilinear Higgs self coupling, because of the absence of the Higgs 
propagator in the contact diagram. 

In this work, we adopt the effective Lagrangian approach, taking the liberty that the 
involved Higgs-boson couplings can be varied freely within reasonable ranges. The relevant 
couplings considered in this work are (i) the top-quark Yukawa coupling, (ii) the trilinear 
Higgs self-coupling, and (iii) the contact-type ttHH coupling. In the top-quark Yukawa and 
contact-type ttHH interactions, we take into account the possibility of the simultaneous 
presence of the scalar and pseudoscalar couplings which can signal CP violation. The ra¬ 
tionale behind the CP-odd part is that the current data, other than the EDM constraints, 
cannot restrict the CP-odd part. The EDM constraints, however, depend on a number 
of assumptions and may therefore be weaken because of cancellation among various CP- 
violating sources im. On the other hand, ttHH is a dimension-5 operator, which may 
originate from a genuine dim -6 operator, e.g., (Qi<I>f/j)(<I>^<h) after electroweak symmetry 
breaking, d) = (0, {v + H)/\/2)^. This operator is thus suppressed only by two powers of 
higher scale, such that it can give a signihcant contribution at the LHC energies. 

Our strategy is hrst to hnd a useful expression for Higgs-boson pair production cross 
sections in terms of these couplings, see Eq. ([^. The coefficient of each term depends 
on the collider energy ^/s, Higgs decay channels, experimental cuts, etc. Thus, such an 
expression enables us to easily obtain the cross section under various experimental conditions 
for arbitrary values of the couplings. Our aim is to extract the information on the Higgs 
couplings, especially on the Higgs self coupling by exploiting the expression. It is helpful 
to consider some experimental cuts which can isolate the contribution from the Higgs self 
coupling and lead to various cross sections with different dependence on the Higgs couplings. 
In this work, specihcally, we employ the 6677 decay mode of the Higgs-boson pair and look 
into the angular separation between the b and b and that between the photons. It is shown 
that one can map out the possible regions of Higgs couplings assuming certain values of 
measured cross sections, though it is channel dependent. Thanks to the largest branching 
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ratio of the Higgs boson into bb, the angular separation between the bottom-quark pair is 
an useful tool in for most of the proposed channel at the LHC. 

In summary, the current work marks a number of improvements over previous published 
works as listed in the following: 

1. We have included the CP-odd part in the top-Yukawa coupling. The CP-even and 
CP-odd parts are constrained by an elliptical-like equation by the current Higgs-boson 
data, as shown in Fig. 4 of Ref. |1]. Note that the effects of the CP-odd part of the 
top-Yukawa coupling at the LHC, ILC, and photon colliders were studied in Ref. IIBI. 
though we study the effects in depth for the LHC here. 

2. We have included the dim-5 anomalous ttHH contact coupling, as done in Ref. [niiii]. 
Furthermore, we also include the CP-odd part of this contact coupling. 

3. We have calculated an easy-to-use expression to obtain the cross sections as a function 
of the involved couplings at each center-of-mass energy. We also obtain similar expres¬ 
sions in various kinematic regions such that one can easily obtain the cross sections 
under the proposed experimental cuts for arbitrary values of the Higgs couplings. 

4. With assumed uncertainties in the measurements of cross sections, we can map out 
the sensitivity regions of parameter space that can be probed at the LHC. 

The organization of the work is as follows. In the next section, we describe the formalism 
for our exploratory approach and present an expression for the Higgs-boson pair production 
cross section in terms of various combinations of the Higgs couplings under consideration. 
In Section HI, we examine the behavior of each term of the cross section versus energies. 
In Section IV, employing the HH —>■ ( 77 ) ( 66 ) decay mode, we illustrate how to extract the 
information on the Higgs couplings by exploiting the angular separations between the Higgs 
decay products. There we also discuss the prospect for the 100 TeV pp machine. We conclude 
in Section V and offer a few comments with regard to our hndings. In Appendix, we compare 
the SM cross sections at 14 TeV with those at 100 TeV for the process pp —)■ iL —)■ 7766 and 
give some comments on the T~^T~bh decay mode. 
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II. FORMALISM 


Higgs-boson pair production via gluon fusion goes through a triangle diagram with a 
Higgs-boson propagator and also through a box diagram with colored particles running in 
it. The relevant couplings involved are top-Yukawa and the Higgs trilinear self coupling. We 
further explore the possibility of a dim-5 anomalous ttHH contact coupling na. They are 
given in this Lagrangian: 


^ (^) + + (1) 

In the SM, = fi'f = 1 and = 0 and = 0. The differential cross section for the 
process g{pi)g{p 2 ) —^ H{p^)H{pi) in the SM was obtained in Ref. [7j as 


d^i.99 HH) 
di 


Gl 




512(27r)3 


X,H9tmFi + {9^tYF, 


S\2 TpSS 
□ 


+ 


i9frGi 


( 2 ) 


where 


D{s) 


3M^ 

s — Fl'jj -|- iMffTjj 


( 3 ) 


and s = (pi P 2 Y, i = (Pi - PsY, and u = {p 2 - PzY pi + p 2 = Ps + P 4 -The loop 
functions = Fa, = Fn, and G^ = Gn with Fa,^ and Gu given in Appendix A.l of 

Ref. [7j. 


Here we extend the result to including the CP-odd top-Yukawa and the anomalous ttHH 
couplings: 


da{gg HH) Gj 


F^s 


dt 


512(27r) 


{XsH9fD{s) + g!Y F| + -f {g^fF, 


pp 

□ 


+ 


{9!YgIY + (pf 


{\.HgYD{^ + gY) Fl + g^YFY,^ ~ + gh^G^^ 


( 4 ) 


More explicitly in terms of each combination of couplings and ignoring the proportionality 


5 






















constant at the beginning of the eqnation, the above eqnation becomes 


ai„ib(s)p [iFina^r + \F^\w] 

+ 2XsH9f^e {D{s) + {F^F^^* + F^Fi^*){g[y] } 

+ + iafr 

+ { 25 Je + |F^"^p + |Gg^p} igfng[r 

+ (gfr 

+ \Fl\\9f,r + \F[\\g2r 

+ 2g2 [XsH9f^e{D)\Fl2 + {gffFlFi^* + {g^fF^F^^*] 

+ 2g2 [A3rf5?e(F)|Flp + F^F^] , (5) 

where F^ = F^, F^^ = Fn, and G^^ = Gn with F^, Fn and Gn given in Appendix A .2 of 
Ref. [ 7 ] while F^^ = Fn, and G^^ = Gn with Fn and Gn in Appendix A .3 of Ref. [ 7 j. In 
the heavy qnark limit, one may have [ 7 j 

Fl = +‘^ + 0 {s/ml), F^^ = -‘^ + 0 {s/ml), F^^ = +‘^ + 0 {s/ml), 

Ff = +1 + 0{s/ml ), F^^^ = -1 + 0{s/ml ), 

G^f = Gg^ = Gg^ = 0{s/ml) ( 6 ) 

leading to large cancellation between the triangle and box diagrams. 

The prodnction cross section normalized to the corresponding SM cross section, with or 
withont cnts, can be parameterized as follows: 

^suigg -t HH) ^ + di{s){g[y] + X3H9t [c2{s){gfy + d2{s){g[)'^] 

+ [c3(s)(^f)^ + d3{s){gff{g2f + di{s){g2Y] 

+X3H [ei{s)g2gft + /i(s)c/f ga] + gft [e2{s){g2Y + f2{s){g2Y] 

+ [e3{s){9ttY + Ms)9fg[ gtt + fYsYg^f] ( 7 ) 

where the nnmerical coefficients 01^2,3(5), <71^2,3,4(5), 61^2,3(5), and /i,2,3,4(s) depend on s and 
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experimental selection cuts. Upon our normalization, the ratio should be equal to 1 when 
gf = = 1 and g[ = gf^^ = 0 or Ci(s) + C 2 (s) + 03 ( 5 ) = 1. The coefficients Ci(s) and 

€ 3 ( 5 ) are for the SM contributions from the triangle and box diagrams, respectively, and the 
coefficient 02 ( 5 ) for the interference between them. 

Once we have the coefficients q, dj, e^, and /ds, the cross sections can be easily obtained for 
any combinations of couplings. Our hrst task is to obtain the dependence of the coefficients 
on the collider energy ^/s, Higgs decay channels, experimental cuts, etc. 


III. BEHAVIOR OF THE CROSS SECTIONS 

We hrst examine the behavior of each piece of cross sections versus energies. We show 
the coefficients c,, dj, e*, /ds at ^/s = 8,14, 33,100 TeV in Table and also in Fig. The 
triangle diagram involves an s-channel Higgs-boson propagator. At center-of-mass energy 
a/I considerably higher than Mh, it behaves like 1/s. Thus the triangle diagram decreases 
more rapidly than the box and contact diagrams, as rehected from the coefficients ci and 
di when ^/s goes from 8 TeV to 100 TeV. The coefficients C3, and associated with 
the box diagram more or less remain the same as the energy goes up. On the other hand, 
the contact diagram involving the coefficients 63 and increases with energy, in particular 
tremendous increases in f^. It is easy to see that the contact diagram is dim-5 and obviously 
grows with energy. At certain high enough energy, it may upset unitarity. 

We can examine the validity of the anomalous ttHH contact coupling by projecting 
out the leading partial-wave coefficient for the scattering ti HH. At high energy, the 
amplitude 

iM{tt ^ HH) ~ 4^ . 

The leading partial-wave coefficient is given by 

= -^ j_^d{cos9)Po{cos9){zM) = • 


^ We note our results are in well accord with those in literature when the comparison is possible. The values 
of Cl, 2,3 and 61 , 2,3 at yT = 100 TeV, for example, are in good agreement with those presented in Ref. M- 
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TABLE I. The behavior of the coefficients Cj(s), di{s), ei{s), and fi{s) versus collider energy ^/s, 
see Eq. Q. 


(TeV) 

Ci(s) 

C 2 (s) 

C3(S) 

di{s) 

\>iH{9t?] 

^2(5) 

\^ZH9l{9ff] 

<* 3 ( 3 ) 

^4(5) 

Mf] 

8 

0.300 

-1.439 

2.139 

0.942 

-6.699 

14.644 

0.733 

14 

0.263 

-1.310 

2.047 

0.820 

-5.961 

13.348 

0.707 

33 

0.232 

-1.193 

1.961 

0.713 

-5.274 

12.126 

0.690 

100 

0.208 

-1.108 

1.900 

0.635 

-4.789 

11.225 

0.683 


(TeV) 

ei(s) 

\^^H9ha] 

62(5) 

\9^{9lf] 

e 3 (s) 

\i9tt?] 

/i(s) 

\hH9[9tt] 

/ 2 (s) 

\9U9ff] 

fsis) 
\ 9 f 9 t 9tt] 

U{s) 

\{9ttf] 

8 

1.460 

-4.313 

2.519 

2.104 

2.350 

-7.761 

3.065 

14 

1.364 

-4.224 

2.617 

1.848 

2.269 

-6.886 

3.769 

33 

1.281 

-4.165 

2.783 

1.622 

2.207 

-6.033 

5.635 

100 

1.214 

-4.137 

2.974 

1.474 

2.154 

-5.342 

10.568 


Requiring |ao| < 1/2 for unitarity we obtain 

Vs < ^ TeV . 

9tt 

Therefore, the anomalous ttHH contact term can be safely applied at the LHC for < 3 — 5 
as most of the collisions occur at \/I < a few TeV. 


To some extent we have understood the behavior of the triangle, box, and contact di¬ 
agrams with the center-of-mass energy, which is kinematically equal to the invariant mass 
Mhh of the Higgs-boson pair. One can then uses Mhh to enhance or reduce the relative 
contributions of triangle or box diagrams. The higher the Mhh the relatively larger pro¬ 
portion comes from the box and contact diagrams. Since Mhh correlates with the boost 
energy of each Higgs boson, a more energetic Higgs boson will decay into a pair of par¬ 
ticles, which have a smaller angular separation between them than a less energetic Higgs 
boson. Therefore, the angular separation ARjj between the decay products is another 
useful kinematic variable to separate the contributions among the triangle, box, and contact 
diagrams. 





























IV. NUMERICAL ANALYSIS 


The Lagrangian in Eq. Q consists of five parameters: the scalar and pseudoscalar parts 

S P 

of the top-Yukawa coupling ’ , the scalar and pseudoscalar parts of the anomalous contact 
coupling , and the Higgs trilinear coupling \^h- In order to facilitate the presentation 
and understanding of the physics, we study a few scenarios: 

1. CPCl-the top-Yukawa coupling involves only the scalar part and the scale in the 
anomalous contact coupling is very large - only g^ and X^h are relevant. The relevant 
coefficients are ci, C 2 , and C 3 . 

2 . CPC2-the top-Yukawa and the anomalous contact couplings involve only the scalar 

part - and Ash are relevant. The relevant coefficients are ci, C 2 , C 3 , ci, 62 , and 

es- 

3. CPVl-the top-Yukawa coupling involves both the scalar and pseudoscalar parts - gf, 
gf, and Ash are relevant. The relevant coefficients are Ci, C 2 , C 3 , di, ^ 2 , d^, and ^ 4 . 

4. CPV2-the contact ttHH coupling involves both the scalar and pseudoscalar parts - 
gfi) 9ttJ Ash are relevant while the top-Yukawa coupling is kept at fixed values. In 
this case, all the coefficients become relevant. In one of the simplest cases with g^ = 1 
and g[ = 0 , for example, the relevant coefficients are Ci, C 2 , C 3 , Ci, 62 , 63 , and f^. 

Note that the above scenarios have been studied using different approaches and separately 
in literature: CPCl in Ref. [12], CPC 2 in Ref. [T3|, and CPVl in Ref. [IB]. The scenario 
CPV2 is new with the CP-odd component of the ttHH, in which the CP-odd component 
gl^ of the ttHH coupling appears either with g^" or in square. 

We used the CTEQ6L1 [19] with both the renormalization and factorization scales /i = 
Mh for the parton distribution functions. Since we focus on the ratios relative to the 
SM predictions, we anticipate the uncertainties due to scale dependence, choice of parton 
distribution functions, and experimental acceptance are reduced to a minimal levelFor the 
branching ratios of the Higgs boson we employ the values for the SM Higgs boson listed in 
the LHC Higgs Cross Section Working Group pni- Here we ignore the slight variation in the 

^ We have observed that the absolute vaues of the production cross sections decrease by about 20 % if we 
take fi = Mhh, well within the theoretical uncertainty estimated in Ref. |S]. 
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diphoton branching ratio due to the change in the top-Yukawa coupling. This is because 
we do not want to interfere with the more important goal of the work - interference effects 
between the triangle and box diagrams and the sensitivity to the trilinear Higgs coupling. 
For Higgs-pair production we used a modihed MADGRAPH implementation [Tnl l23] , which 
allows us to vary the top-Yukawa, Higgs trilinear, and the contact ttHH couplings. 

In this work, we make a working assumption that the 6677 background can be estimated 
with a reasonable accuracy and can be extracted from the experimental data. Once the 
background is subtracted from the data, we are left with the signal events. In order to do 
so, we impose a set of basic cuts for the acceptance of the b quarks and photons, and use 
R-tagging according to the ATLAS template in Delphes v.3 m, as well as an invariant 
mass window on the bb pair and the photon pair around the Higgs boson mass. The basic 
cuts are 

Pt^ > 10 GeV, pt^ > 10 GeV, \p^\ < 2.5, \pb\ < 2.5, > 0.4, AR^b > 0.4, 

< 15GeV, |Mh - < 25 GeV, |Mh - < 25 GeV. (8) 

Note that in the appendix we shall discuss the feasibility of using the T~^T~bh mode. With 
this set of basic cuts we continue to study the signal events in various kinematical regions 
separated by AR^^ and ARbb- 

A. CPCl: gf and X3H 

This is the simplest scenario to investigate the variation of the triangle and box diagrams 
with respect to changes in the Yukawa and trilinear couplings. The corresponding coefficients 
in Eq. ([^ are Ci, C 2 , C 3 . We let the Higgs boson pair decay into 

HH ( 77 ) (bb) . 

Note that if the process is studied without detector simulation, the distributions for other 

decay channels, like bhrT or yyrr would be the same. Nevertheless, the resolutions for b, r, 

^ According to our most recent analysis using all the Run I data [4], the allowed range of the top-Yukawa 
coupling is, at 95% C.L., [—1.1,—0.8] U [0.7,1.4] in the scenario where gauge-Higgs coupling and top-, 
bottom-, and tau-Yukawa couplings are varied; and additionally if we included non-SM contributions to 
Hgg and Tfyy vertices, the range becomes [—2, 2]. The values of top-Yukawa coupling used in this study 
are in the allowed ranges such that the diphoton signal strength is kept within the uncertainty of the 
experimental data, for example, HggH = 1-32 ± (ATLAS) [21] and HggH = 1-12 io .32 (CMS) [22] for 
gg^ H ^ 77 . 



and 7 are quite different in detectors, and so are the backgrounds considered for each decay 
channel. In the following, we focus on the 7766 channel, which has been considered in a 
number of works. 

We use MADGRAPH v .5 [ 21 ] with parton showering by Pythia v.6 [ 21 ], detector simula¬ 
tions using Delphes v. 3 [ 2 S], and the analysis tools by MadAnalysis 5 [ 2 Z|. We have verihed 
that the coefficients ci,2,3 using the default, ATLAS, and CMS templates inside Delphes v .3 
are within 10 % of one another. From now on we employ the ATLAS template in detector 
simulations. We show the coefficients in Table |TT 1 We will come back to this table a bit later. 

We hrst look at the distribution of events versus the invariant mass of the Higgs-boson 
pair via HH ^ 7766 in Fig. (upper panel) and versus the transverse momentum of the 
photon pair (lower panel). In principle, up to detector simulations and higher-order 
emissions, the transverse momentum distribution of the photon pair is the same as that of 
the bb pair. The triangle diagram (red line) peaks at the lower invariant mass and decreases 
with because of the s-channel Higgs propagator. The box diagram (skyblue line), on 

the other hand, is larger than the triangle diagram at high invariant mass. The (darkblue) 
SM line represents the whole contributions including the destructive interference between 
the triangle and box diagrams. Similar behavior can be seen in the pt distribution of the 
photon pair. 

Next, we look at the distributions versus and ARhh in Fig. The lines are the 

same as in Fig. [^ As we have explained in the previous section, AR^^ and AR^b between 
the decay products of each Higgs boson are useful variables to separate the triangle and 
box contributions. The angular distribution AR between the two decay products of each 
Higgs boson correlates with the energy of the Higgs boson, which in turns correlates with 
the invariant of the Higgs-boson pair. The higher the invariant mass, the more energetic 
the Higgs boson will be, and the smaller the angular separation between the decay products 
will be. Therefore, the triangle diagram has wider separation than the box diagram. 

It is clear that the distributions of AR^^ and ARbb have similar behavior within uncer¬ 
tainties. The box diagram and also the SM, which is dominated by the box contribution, 
have a peak at ARj^ or ARbb less than 2.0, while the triangle diagram prefers to have the 
majority at larger AR^^ or ARbb, say between 2 and 3 . We therefore come up with (i) 
ARyy > 2 (< 2), (ii) ARbb > 2 (< 2), and (hi) AR^.y > 2 and ARbb > 2 (both < 2) to 
enrich the sample of triangle (box) contribution. In the following, we use AR to denote 
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either or ARbb, unless stated distinctively. 

We can now look at Table where the coefficients for the ratio of the cross sections 
a{gg -)■ HH)/asM{gg ^ HH) as in Eq. 0 are shown. In the CPCl case, the relevant 
coefficients are ci, C 2 , and C 3 in which ci is induced by the triangle diagram, C 3 by the box 
diagram, and C 2 by the interference between them. The rows labeled “Basic Cuts” are the 
ratio of cross sections under the set of cuts in Eq. (|^. In the same Table, we also show 
the coefficients obtained after applying the angular-separation cuts of AR > 2 or < 2 and 
both ARyy^bb > 2 or < 2 . It is clear that AR-y-y > 2 (< 2 ) enriches the triangle-diagram 
(box-diagram) contribution. Similar is true for AR^b > 2 (< 2). Further enhancement of 
triangle diagram can be obtained with both AR^^ > 2 and ARbb > 2 , and vice versa for box 
diagram. 

In the following, we investigate the sensitivity in the parameter space {g^, Xsh) that 
one can reach at the 14 TeV with 3000 fb“^ luminosity by using the measurements of cross 
sections in various kinematical regions. Since we have found that the triangle and box 
contributions can be distinguished using the AR cuts, we make use of the measured cross 
sections in the kinematical regions separated by these cuts. 

There are two issues that we have to considered when we take the measured cross sections 
in the kinematical regions. First, the SM backgrounds for the decay channel that we consider, 
and second, the Next-to-Leading-Order (NLO) corrections [28]. It was shown in Ref. [28| 
that the NLO and NNLO corrections can be as large as 100% with uncertainty of order 
10-20%. The SM backgrounds, on the other hand, can be estimated with uncertainties less 
than the NNLO corrections. We therefore adopt an approach that the signal cross sections 
(after background subtraction) are measured with uncertainties of order 25-50%. 


About the signal cross sections, the hrst and the second columns of Table |111| show 

the SM cross sections for the process pp — HH —)■ 7766 with detector simulations under 

various cuts at the 14 TeV LHC 0 We have taken account of the SM NLO cross section 

o'suipp —t HH) ~ 34 fb, the Higgs branching fractions, and both the photon and 5-quark 

reconstruction efficiencies with angular-separation cuts of ARbb,-yy In various kinematical 

^ Before applying the basic cuts, we find the cross sections are 8.92 x 10“^ and 2.41 in fb for SM-14(7766) 
and SM-14(T+r“66), respectively, which agree very well with those in Ref. [S]. While our cross sections 
after applying the basic cuts are smaller than those presented in Ref. [8| by a factor of ~ 3 for 77 && and 
a factor of ^ 30 for This is basically because we have implemented full detector simulation to 

reconstruct b quarks, photons, and r leptons from Higgs and partly due to different experimental cuts 
applied and different b- and r-tagging efficiencies taken. One may need to optimize the cuts to increase 
the signal to background ratio but it is beyond the scope of this paper and we will pursue this issue later 
in our future publication. Incidentally the SM-lofe^yyfefe) cross section is 3.73 fb before applying the basic 
cuts. 




regions depending on the angular cuts, the cross sections range from ~ 0.001 fb to ~ 0.01 
fb. With an integrated luminosity of 3000 fb~^, we expect of order 30 signal events when the 
cross section is 0.01 fb. An estimate of the statistical error is given by the square root of the 
number of events \/iV, which is then roughly 20% of the total number. Taking into account 
the uncertainty of order 10 — 20% from NLO and NNLO corrections, in this work, we use 
a total uncertainty of 25 — 50% in the signal cross section in the estimation of sensitivity 
of the couplings. Our approach is more or less valid except for the case in which both 
the ARhh > 2 and > 2 cuts are imposed simultaneously. It would be challenging 

to measure this size of cross section only in the HH ^ 7766 mode and one may need to 
combine the measurements in different Higgs-decay channels. Or, one may rely on the future 
colliders such as a 100 TeV pp machine with larger cross sections and/or higher luminosities. 

In Fig. 1^ we show the contour lines of a{gg — HH)/asM{gg —t HH) = 0.5 and 1.5 with 
the Higgs boson pair decaying into ( 77 ) ( 66 ). In each panel, we assume three measurements 
of the ratios corresponding to basic cuts (orange lines), AR > 2 (dashed black lines), and 
AR < 2 (solid black lines): here AR represents AR^^ (upper left), AR^i, (upper right), or 
AR.y^^bb (lower). Therefore, for example, if the bacis-cuts cross section ratio is measured to 
be consistent with the SM prediction within 50% error, any points in the two bands bounded 
by the two pairs of orange lines are allowed. In each band, a rather wide range of g^ and 
is allowed although they are correlated. Suppose we only make one measurement of 
the cross section without or with a cut on AR, we would not be able to pin down useful 
values for g^ and X^h- However, since the shape of the three bands are not exactly the 
same, we can make use of three simultaneous measurements in order to obtain more useful 
information for the couplings gf and X^h- 

In the upper-left panel of Fig. we suppose that one can make three measurements of 
cross sections: with basic cuts, AR^^ > 2, and AR^^ < 2. We assume that the measurements 
agree with the SM predictions within 25% or 50% uncertainty. The region of parameter space 
in [gf, X^h) bounded by the three measurements is shown by the lighter purple region for 
50% uncertainty and darker purple region for 25% uncertainty. Similarly the upper-right 
panel is for the regions with the ARbt cut. In the lower panel, we show the regions with 
the combined cuts of AR^^ and ARbb'- both larger than or smaller than 2. The implications 
from the measurements are very significant. First, all panels show that gf is significantly 
away from zero if one can simultaneously measure the cross sections (no matter with 25% 
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or 50% uncertainties) with basic cuts, > 2, and < 2; and similarly for ARi,b 

and using both distributions. Second, as shown in the lower panel, the value for is 
statistically distinct from zero if one measures the cross sections with a 25% uncertainty. 
This is achieved by using both AR^^ and AR^b > 2 or < 2. Furthermore, from the lower 
panel in Fig. we can see that with 25% level uncertainty the values of sensitivity 
regions are 0.3 < lAs^l < 2.6. 

We can repeat the exercise with the measured cross sections being multiples of the SM 
predictions. We show the corresponding 25% and 50% regions in Fig. I for cr/cTsM = 
0.5, 1, 2, 5, 10. Only with both AR^^ and ARbb, one can really tell if X^h is signihcantly 
distinct from zero. The sensitivity regions for X^h for various ct/ctsm are indicated by the 
darker color areas. 


If the top-Yukawa coupling can be constrained more effectively by Higgs production, 
by tin production, or by single top with Higgs production in the future measurements, 
say gf = 1 ± 0.1 (10% uncertainty), it can help pinning down the acceptable range of 
X^H- However, even in this case, we emphasize the importance of simultaneous independent 
measurements, as illustrated in the following argument. In the limit of = 1, the ratio of 
the cross sections is given by 


ajgg ^ HH) 
(^suigg -t HH) 


— Cl(s) Ag^ + C2(s) XsH + ^ 3 ( 3 ) 


(9) 


Suppose (j{gg HH) is measured to be the same as crsuigg —t HH) and then, using the 
relation Ci(s)+C2(s)+C3(s) = 1, one may hnd the two solutions for As/;-: 1 or — C2(s)/ci(s) —1. 
For example, one may have X^h = 1 or 4 at most if only the basic-cuts ratio is measured, see 
Table [Hj Therefore, one cannot determine X^h uniquely with only one measurement even 
when the measurement is very precise and the exact value of gf is known. It is unlikely to 
resolve this two-fold ambiguity at the LHC even we assume the three measurements of the 
ratios, as shown in Fig. Also, the situation remains the same at the 100 TeV pp machine 
in which we have X^h = 1 or 5 in the bacis-cuts case when a{gg — HH) = asuigg — t HH), 


see Table IV and Fig. 24 If a future e'''e linear collider and/or the 100 TeV pp machine are 


operating in the era of the high-luminosity LHC, combined efforts are desirable to determine 
the value of X^h uniquely 


14 





B. CPC2: g^, X^h, and 


This is the scenario that involves all scalar-type couplings in the triangle, box, and contact 
diagrams. The corresponding coefficients in Eq. ([^ are 01 , 02 , 03 , 61 , 62 , 63 . Results at the 
detector level using the ATLAS template in Delphes v.3 are shown in Table |TTj 

We hrst examine the cross section versus one input parameter at a time, shown in Fig. 
while keeping the two parameters at their corresponding SM values. In the upper-left panel 
for a/asM versus Xsh, the lowest point occurs at X^h ~ 2.5 when the interference term 
strongly cancels the triangle and box diagrams. Then the ratio increases from the lowest 
point on either side of X^h ~ 2.5. Negative X^hs give constructive interference while positive 
X^hs give destructive interference. One may observe similar behavior when is varied as 
shown in the lower panel. Taking X^h = S'f = 1, 

Since 61 (s) > 0 and 62 ( 5 ) < 0, we see that the contact diagram interferes constructively with 
the triangle diagram but destructively with the box diagram. The dominance of the box 
diagram leads to the totally destructive interference when > 0 , resulting in the minimum 
at g^t ^ 0.5. 

We show contours for the ratio cr/asM = 0.5,1.5 in the plane of (A 3 i^, g^) (upper-left), 
i^ 3 H, 9u) (upper-right), and {gf, gf^) (lower) in Fig. 0 The dashed lines denoted by —50% 
is for (t/ctsm = 0.5 and those by -|-50% for a/asu = 1-5. In the upper-left panel in the 
plane of (X^h, gt)i we show contours for gf^ = 0,1. The ( 7 ^ = 0 is the same as the SM so 
that the contours are exactly the same as in Fig. while with ( 7 ^ = 1, the contact diagram 
contributes signihcantly to the cross section, so that the contours shift more to negative 
(positive) X^H for positive (negative) g^. In the upper-right panel, where we show the 
contours in the plane of (A 3 //, 77 ^^), the gf^ negatively correlates with X^h because 61 (s) > 0 . 
In the lower panel, where we fix X^h = 1, somewhat nontrivial correlation between g^ and 

exists. 

We use the same tools as in the CPCl case to investigate the decay channel HH —7766 
with parton showering and detector simulations. We hrst show the invariant mass 
distribution of the Higgs-boson pair and the transverse momentum distribution of the photon 
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pair in Fig. The figure clearly illustrates the behavior of each contributing diagram. The 
triangle diagram (red lines) peaks at the lower invariant mass and decreases with 
because of the ,s-cha nn el Higgs propagator. Then followed by the box diagram (skyblue 
lines) which is larger than the triangle diagram at high invariant mass. The (darkblue) SM 
line represents the contributions from the triangle and box diagram including the destructive 
interference between them. These three distributions are the same as in the CPCl case. The 
contact diagram (green lines) shows similar behavior as the box diagram (skyblue lines) at 
low and but with higher and larger tails. The grey and magenta lines represent 

the full contributions including the destructive and constructive interferences among the 
triangle, box, and contact diagrams when = +1 and gf^ = —1, respectively. The contact 
diagram is the largest at the high invariant mass. It demonstrates what we describe earlier 
that the contact diagram grows with energy. 

We show the angular distributions and ARbh between the two decay products of 

each Higgs boson in Fig. The lines are the same as in Fig. Similar to the CPCl 
case, the higher the invariant mass, the more energetic the Higgs boson, and the smaller 
the angular separation between the decay products will be. Therefore, triangle diagram 
(red lines) has the widest separation, then followed by the box diagram (skyblue lines), and 
finally the contact diagram (green lines) has the smallest angular separation. We come up 
with the similar cuts as in the CPCl case: AR larger or smaller than 2 to discriminate the 
triangle, box, and contact diagrams. We show in Table |IT| the coefficients Ci, C 2 , C 3 , Ci, 62 , 63 
such that the ratio of cross sections to the SM predictions can be given by Eq. ([^. 

Similar to what we have done for CPCl, we can make use of three simultaneous mea¬ 
surements of cross sections with basic cuts, AR > 2, and AR < 2. We show the region of 
parameter space that we can obtain using AR^^ (upper panels), AR^b (middle panels), and 


AR^^ and ARbb (lower panels) in the plane of 9t) Fig- 10 Those on the left are for 
g^^ = 1 while those on the right are for g^^ = —1. Similarly, we show the parameter space in 


the plane of {X^Hi 9t) i^^ Fig. 11 and in the plane of {g^, g^^) in Fig. 12 


C. CPVl: gf, gP, and X^ 

In this scenario, we entertain the possibility that the top Yukawa coupling allows an 
imaginary part. In most of the measurements of the Higgs boson production cross sections. 
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for example, Higgs boson production cross section via gluon fusion and ttH production, both 
the real and imaginary parts of the coupling come in the form \gf\‘^ + therefore one 

cannot tell the phase in the coupling]^ The relevant coefficients for this CPVI scenario are 
Cl, C 2 , C 3 , di, ^ 2 , ds, di- They are shown in Table [ h] at the detector simulation level (ATLAS). 

We first show the variations of cross sections versus X^h with some fixed values of and 
gf in Fig. m Also, the contours for the ratio ct/ctsm = 1 in the plane of {Xsh, 5'f) (upper- 
left), {X 3 H, gf) (upper-right), and {g^, gf) (lower) for a few values of the third parameter 
are shown in Fig. 

Similar to previous two scenarios, we use the same tools to analyze the decay channel 
HH —)■ 7766 with parton showering and detector simulations. We show the invariant mass 


and in Fig. 15, and the angular distributions and Ad?bf, between the two 

decay products of each Higgs boson in Fig. [T^ 

The terms by the triangle diagram (proportional to Ci and di in red and orange lines, 
respectively) give the widest separation among all the terms. The terms by the box diagram 
(proportional to C 3 and in skyblue and blue lines, respectively) give smaller angular 
separation. The full set of diagrams at the SM values (darkblue lines) and at gt = gf = 1/ 
(grey lines) give similar results as the box diagram. 

Similar to the CPCl and CPC2 cases, we use the cuts AR larger or smaller than 
2 to discriminate the triangle and box diagrams. We show in Table the coefficients 
Cl, C 2 , C 3 , di, d 2 , ^ 3 ,(^ 4 , which are relevant ones in the CPVl scenario, such that the ratio 
of cross sections to the SM predictions can be obtained by Eq. (|^. We show the region of 
parameter space that we can obtain using AR^^, ARbt, and AR^^ and ARbt in the plane 


of (^ff, gf) in Fig. 17, in the plane of {Xsh, gf) Fig. 18, and in the plane of {X 3 H 1 gt) 
Fig. Ill 


D. CPV2: and X^h 

Here we study another CP-violating scenario with the CP-even and CP-odd components 

of the ttHH coupling with the top-Yukawa couplings g^ and gf at fixed values. Note that 

the CP-odd coupling gj^ only appears in the product with gj^ or by itself squared. In this 

® Single top plus Higgs production, on the other hand, has some chances to isolate the phase of the Yukawa 
coupling [50] . 

® We note that the lines are the same for the negative values of gf since the cross section contains only the 
(gf )2 and (gf )^ terms. 
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case, all the coefficients are relevant and they are shown in Table [IT] at the detector simulation 
level (ATLAS). 


Similar to previous scenarios we used HH ^ 7766 with parton showering and detector 
simulations. We use the cuts AR larger or smaller than 2 to discriminate the triangle and 
box diagrams. Using the coefficients presented in Table the ratio of cross sections to the 
SM predictions can be obtained by Eq. ([^. We show the region of parameter space that we 
can obtain using AR^^, AR^b, and AR^^^bb in the plane of (( 7 ^, for hxed X^h = 1 , fi'f = 1 , 

and gf = 0 in Fig. 20 and similarly in the plane of {Xsh, gtt) 9t — 9t = 9tt — ^ 
Fig. 21 and hnally in the plane of (Ash, g^^) for g^ = 1, g^ = 0, and g^ = 0.5 in Fig. 
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E. 100 TeV Prospect 


All the results represented for the 14 TeV run were obtained by manipulating the coef- 
hcients represented in Table We represent the coefficients Ci_ 2 , 3 , ei, 2 , 3 ) /i, 2 , 3,4 for 


the 100 TeV pp machine in Table IV Just for illustrations, we show the distributions of 
the invariant mass M^^bl ^^nd angular separation AR^^ for the CPCl case at the 100 TeV 


machine in Fig. ^ We found that the behavior of the distributions at 100 TeV is very 
similar to those at 14 TeV. Therefore, the kinematic regions of interests separated by AR 
can be taken to be the same as 14 TeV. We can make simultaneous measurements of cross 
sections at 100 TeV pp machine to isolate the Higgs trilinear coupling. We show the sensi¬ 


tivity regions of parameter space in the CPCl case at the 100 TeV pp machine in Fig. 24 


The regions are very similar to those in 14 TeV, though not exactly the same. Sensitivity 
reach for each coupling in other cases can be obtained by similar methods with the assumed 
luminosity. 


V. CONCLUSIONS 

In this work, we have studied the behavior of Higgs-boson pair production via gluon fusion 
at the 14 TeV LHC and the 100 TeV pp machine. We have performed an exploratory study 
with heavy degrees of freedom being integrated out and resulting in possible modihcations of 
the top-Yukawa coupling, Higgs trilinear coupling, and a new contact ttHH coupling, as well 
as the potential CP-odd component in the Yukawa and contact couplings. We have identihed 
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useful variables - the angular separation between the decay products of the Higgs boson - 
to discriminate among the contributions from the triangle, box, and contact diagrams. We 
have successfully demonstrated that with three simultaneous measurements of the Higgs- 
pair production cross sections, dehned by the kinematic cuts, one can statistically show a 
nonzero value for the Higgs trilinear coupling X^h if we can measure the cross sections with 
less than 25% uncertainty. This is the key result of this work. 

We also offer the following comments with regards to our hndings. 

1. The triangle diagram, which contains an s-channel Higgs propagator, does not increase 
as much as the box diagram or the contact diagram with the center-of-mass energy 
y/S. This explains why the opening angle {AR^^ or ARbb) in the decay products of 
each Higgs boson is a useful variable to separate between the triangle and the box 
diagram. Thus, it helps to isolate the Higgs trilinear coupling X^h- 

2. The contact diagram contains a dim-5 operator ttHH, which actually breaks the 
unitarity at about y/S 17.6/(71 TeV. This implies that it could become dominant at 
high invariant mass. 

3. Suppose we take a measurement of cross sections, we can map out the possible region 
of parameter space. Since in different kinematic regions the regions of parameter space 
are mapped out differently, such that simultaneous measurements can map out the in¬ 
tersected regions. With measurement uncertainties less than 25% one can statistically 
show a nonzero value for the Higgs trilinear coupling, and also obtain the sensitivity 
regions of Xsh- 0.3 < lAs^l < 2.6. for a/asu = 1- 

4. We found that the behavior of the distributions of the invariant mass and angular 
separation AR^^ or ARb^ at 14 TeV are very similar to those at 100 TeV. We can then 
use the same method as in 14 TeV to isolate the Higgs trilinear coupling. 

5. It is difficult, if not impossible, to determine the Higgs trilinear coupling uniquely 
at the LHC and 100 TeV pp machine even in the simplest case assuming very high 
luminosity and precise independent input for the top-Yukawa coupling. We suggest 
to combine the LHC results with information which can be obtained at a future 
linear collider. 
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6 . If the couplings deviate from their SM values, the Higgs-boson pair production cross 
section can easily increase by an order of magnitude. For example, in the CPC2 case, 
c’‘/o‘SM > 10 for > 9or < —4 when gf = 1 and = 0, gf > 1.7or < —1.3 when 
X^H = 1 and g^-i- = 0, and gf^ > 2.6 or < —1.4 when X^h = Qt = 1^ see Fig. The 
cross section larger than the SM prediction may reveal the new physics hidden behind 
the SM and we can have better prospect to measure the Higgs self coupling at the 
LHC. 
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Appendix A 


In Table 111, we show the SM cross sections for pp —HH —7766 at the 14 TeV LHC 
with and without angular-separation cuts. Note that the cross section before applying any 
cuts is about 0.09 fb and it becomes 0.005 fb after applying the basic cuts. In the region 
of > 2(< 2), the cross section is 0.0013 fb (0.0038 fb) where it is dominated by the 

triangle (box) diagram. The ratio is about 1 : 2.8. We also show the cross sections for the 
100 TeV pp machine, and the corresponding ratio is about 1 : 3.7. It shows the fact that the 
triangle diagram is more suppressed because of the s-channel Higgs propagator at higher 
energy. In the regions of ARhh larger and smaller than 2, the ratios are 1 : 5.7 and 1 : 8.2 at 
the 14 TeV LHC and the 100 TeV pp machine, respectively. 

As we have promised, we are going to comment on the HH — T^T~bb decay mode. This 
mode has the obvious advantage of a larger branching ratio than the 7765 mode, but the 
identihcation efficiency and momentum measurements of r leptons are much weaker than 


photons. In Table we show the SM cross sections for pp —)■ HH —r+r bb at the 14 
TeV LHC with and without angular-separation cuts. Taking into account the branching 
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ratios, and the identification and selection efficiencies, the event rates of bb are similar 
to those of the 'yybb mode. 
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TABLE II. 14 TeV LHC: The coeffcients for the ratio of the cross sections (T{gg —>■ 
HH)/(TsM{gg HH) as in Eq. 0 with and without the angular-separation cuts of > 

or < 2; ARhh > or < 2; and AR^^ and ARhh both > 2 or < 2. The relevant coefficients for the 
CPCl scenario are ci, C 2 , cs; those for the CPC2 scenario are ci, C 2 , 03 , 61 , 62 , 63 ; and those for 
the CPVI scenario are ci, 62 , C 3 ,di,d 2 ,d 3 ,di. All the coefficients are involved in the CPV2 case. 
Results are at the detector level using the ATLAS template in Delphes v.3. 


^ : 14 TeV 


ci(s) 

C 2 (s) 

C 3 (s) 

di{s) 

d2{s) 

dsis) 

d^i^s) 

Cuts 



1 tjf)-! h 

Inigm [>^3Hgf{gm Mfigry 


No cuts 


0.263 

-1.310 

2.047 

0.820 

-5.961 

13.348 

0.707 

Basic Cuts 


0.221 

-1.104 

1.883 

0.665 

-4.738 

11.757 

0.650 

AR^^ > 2 


0.470 

-1.868 

2.398 

1.481 

-9.754 

19.859 

0.858 

< 2 


0.133 

-0.834 

1.701 

0.376 

-2.959 

8.884 

0.576 

ARbb > 2 


0.666 

-2.512 

2.847 

2.040 

-13.425 

25.316 

1.074 

ARbb < 2 


0.143 

-0.857 

1.714 

0.424 

-3.214 

9.378 

0.575 

ARbb > 2 & ARryry 

> 2 

0.895 

-3.150 

3.255 

2.613 

-17.210 

30.456 

1.278 

ARbb < 2 & AR^^ 

< 2 

0.121 

-0.785 

1.664 

0.319 

-2.630 

8.257 

0.563 

: 14TeV 


ei(s) 

62(5) 

63(5) 

/i(s) 

/ 2 (S) 

/3(s) 

his) 

Cuts 


\^3Hgfgtt\ 

Igftigff 

1 M?] [A 3 ^.«] 

\gttigff] 

gft] 

liguf] 

No cuts 


1.364 

-4.224 

2.617 

1.848 

2.269 

-6.886 

3.769 

Basic Cuts 


1.381 

-3.966 

2.521 

1.939 

2.328 

-5.239 

3.178 

AR^^ > 2 


1.857 

-4.506 

2.267 

4.014 

2.555 

-11.188 

2.569 

AR^^ < 2 


1.212 

-3.774 

2.611 

1.203 

2.247 

-3.130 

3.394 

ARbb > 2 


2.248 

-5.214 

2.474 

5.517 

3.367 

-16.349 

3.003 

ARbb < 2 


1.229 

-3.747 

2.529 

1.311 

2.146 

-3.290 

3.208 

ARbb > 2 & AR^^ 

> 2 

3.047 

-5.947 

2.780 

7.274 

3.759 

-21.142 

3.547 

ARbb < 2 & AR^^ 

< 2 

1.238 

-3.758 

2.664 

1.095 

2.211 

-2.716 

3.500 
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TABLE III. The SM cross sections for the process pp ^ HH ^ 7766 with various angular- 
separation cuts on at the 14 TeV LHC (second column) and at the 100 TeV pp machine 

(third column). The last column shows them for the process pp —5- HH — T^T~bb at the 14 
TeV LHC. We have taken account of the SM NLO cross section a^uipp —)• LlLl) ^ 34 fb, the 
Higgs branching fractions, and both the photon and 6 -quark reconstruction efficiencies. The pT 
dependence of 6 -tagging efficiency is considered and 0.5 is taken for the r-tagging efficiency. Also 
considered is the mis-tagging probability of Pj^r = 0.01. Results are at the detector level using 
the ATLAS template in Delphes v.3. 



Cuts 


SM-14 ( 7766 ) 

SM-100 ( 7766 ) SM-14 

Cross Section (fb) 

(p 

*-r bb) 


No cuts 


8.92 

X 

10 - 

-2 

3.73 




2.41 





Basic Cuts 


5.1 X 10“=^ 


2.05 

X 

10 - 

-1 

3.53 

X 

10 - 

-3 


^ 2 


1.34 

X 

10 - 

-3 

4.34 

X 

10 - 

-2 

6.43 

X 

10 - 

-4 


j 2 


3.76 

X 

10 - 

-3 

1.61 

X 

10 - 

-1 

2.89 

X 

10 - 

-3 


ARbb > 2 


7.61 

X 

10 - 

-4 

2.23 

X 

10 - 

-2 

4.82 

X 

10 - 

-4 


ARbb < 2 


4.34 

X 

10 - 

-3 

1.83 

X 

10 - 

-1 

3.05 

X 

10 - 

-3 

ARbb 

> 2 & AR^^iT.+r- 

> 2 

4.79 

X 

10 - 

-4 

1.40 

X 

10 - 

-2 

3.21 

X 

10 - 

-4 

ARbb 

< 2 & AR^^It+^- 

< 2 

3.48 

X 

10 - 

-3 

1.53 

X 

10 - 

-1 

2.73 

X 

10 - 

-3 
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TABLE IV. 100 TeV pp machine: The coeffcients for the ratio of the cross sections cr{gg —>■ 
HH)/asM{gg HH) as in Eq. 0 with and without the angular-separation cuts of > 

or < 2; ARhh > or < 2; and AR^^ and ARhh both > 2 or < 2. The relevant coefficients for the 
CPCl scenario are ci, C 2 , cs; those for the CPC2 scenario are ci, C 2 , 03 , 61 , 62 , 63 ; and those for 
the CPVI scenario are ci, 62 , C 3 ,di,d 2 ,d 3 ,di. All the coefficients are involved in the CPV2 case. 
Results are at the detector level using the ATLAS template in Delphes v.3. 



\/~s : 100 TeV 



ci(s) 

C2(s) 

C3(s) 

di{s) 

d2is) 

dsis) 

dhs) 


Cuts 



CO to 

to 

IhnigfY 

1 tjf)-! h 

InigHY 

\^3Hgfig[)Y 


1 f(9,n‘i 


No cuts 



0.208 

-1.108 

1.900 

0.635 

-4.789 

11.225 

0.683 


Basic Cuts 



0.173 

-1.032 

1.860 

0.503 

-4.045 

10.019 

0.633 


AR^^ > 2 



0.389 

-1.904 

2.515 

1.275 

-6.972 

13.375 

0.853 


AR^^ < 2 



0.115 

-0.798 

1.683 

0.295 

-3.258 

9.116 

0.574 


ARbb > 2 



0.607 

-2.419 

2.813 

1.845 

-9.336 

17.393 

1.057 


ARbb < 2 



0.120 

-0.863 

1.743 

0.340 

-3.400 

9.119 

0.581 

Ai?, 

> 2 & ARryry 

> 

2 

0.753 

-2.662 

2.909 

2.248 

-10.518 

17.691 

1.245 

Ai?, 

!,{, < 2 & ARryry 

< 

2 

0.102 

-0.733 

1.632 

0.249 

-3.041 

8.700 

0.565 


^/s : lOOTeV 



ei(s) 

62 ( 5 ) 

63 ( 5 ) 

his) 

his) 

his) 

his) 


Cuts 



\^3Hgfgtt\ 

IgftigfY 

1 MY] [^sHg^gYb] 

IguigHY 

\gtgY gu] 

\igtt)Y 


No cuts 



1.214 

-4.137 

2.974 

1.474 

2.154 

-5.342 

10.568 


Basic Cuts 



1.170 

-4.081 

2.848 

1.300 

1.935 

-3.379 

7.802 


AR^^ > 2 



1.782 

-4.886 

2.591 

3.675 

2.151 

-2.696 

5.511 


AR^^ < 2 



1.006 

-3.865 

2.917 

0.662 

1.878 

-3.563 

8.419 


ARbb > 2 



2.011 

-5.585 

2.957 

6.947 

2.576 

-4.961 

5.373 


ARbb < 2 



1.068 

-3.898 

2.834 

0.612 

1.857 

-3.186 

8.099 

AR, 

> 2 & AR^^ 

> 

2 

2.483 

-5.858 

3.106 

8.165 

2.694 

-4.722 

6.079 

AR, 

< 2 & AR^^ 

< 

2 

0.995 

-3.798 

2.928 

0.437 

1.851 

-3.466 

8.647 
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FIG. 1. The coeffcients Cj(s) and ej(s) (upper) with i = 1,2,3 and dj{s) and fj{s) (lower) with 
j = 1, 2, 3,4 as functions of s. 
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LHC-14, Detector Level-ATLAS 



LHC-14, Detector Level-ATLAS 



FIG. 2. CPCl: Distributions in (upper) and (lower) in the decay products of the 

Higgs boson pair with detector simulation. 















LHC-14, Detector Level-ATLAS 



LHC-14, Detector Level-ATLAS 



FIG. 3. CPCl: Distributions in AR^-y (upper) and ARhi, (lower) in the decay products of the 
Higgs boson pair with detector simulation. 
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FIG. 4. CPCl: The 25% and 50% sensitivity regions bounded by three measurements of 
cross sections with basic cuts, > 2, and AR^^ < 2 (the upper-left panel); with basic 

cuts, > 2, and < 2 (the upper-right panel); and basic cuts, AR^^, ARhi, > 2, and 

AR^ry, AiZfcfe < 2 (the lower panel). We assume that the measurements agree with the SM values 
with uncertainties of 25% and 50%, respectively. 
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FIG. 5. CPCl: The same as Fig. |^but we assume ct/ctsm = 0.5 — 10. 
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FIG. 6. CPC2: Variation of the ratio a/a^M using Eq. Q versus \^h (upper-left), (upper- 
right), and gfi (lower) while keeping the other two parameters at their corresponding SM values. 
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g", 8"« 




- 3 - 2-10 1 2 3 


FIG. 7. CPC2: Contours marked by (t/itsm = 0.5,1.5 for {X^h, 9t) (upper-left), {X^h, 9tt) 
(upper-right), and (gf, gf^) (lower). 
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LHC-14, Detector Level-ATLAS, CPC, ?i3j^=g®=1 



LHC-14, Detector Level-ATLAS, CPC, >.3H=g®=1 



FIG. 8. CPC2: Distributions in (upper) and (lower) in the decay products of the 

Higgs boson pair with detector simulation. 
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LHC-14, Detector Level-ATLAS, CPC, 



LHC-14, Detector Level-ATLAS, CPC, >.3„=g^^=1 



FIG. 9. CPC2: Angular distributions of and AR^i, between the two photons and between 

the two b quarks. 


35 













































25% & 50% lortr/cTsM = 1 

gstt = -1 

Detector Level - ATLAS, LHC-14 

Basic Cuts - 

ARhb >2 - - - 

ARbb < 2 




ss- 


25% & 50% f'ortr/0-sM = 1 
gstt = -1 

Detector Level - ATLAS, LHC-14 

Basic Cuts - 

AR^y^bb ^ 2 
^J^j')'.bb < 2 


SP' 

L—.—^—■ -‘i ■—■—■—^— 





FIG. 10. CPC2: The 25% and 50% sensitivity regions in {XsH^dt) bounded by three mea¬ 
surements of cross sections with basic cuts, AR.yy > 2, and < 2 (the upper panels); with 

basic cuts, > 2, and AR^i, < 2 (the middle panels); and basic cuts, AR^y^, ARhb > 2, and 

ARy^, ARbb < 2 (the lower panels). The left panels are for = 1 while those on the right are for 
= —1. We assume that the measurements agree with the SM values with uncertainties of 25% 
and 50%, respectively. 
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FIG. 11. CPC2: similar to Fig. 


but in the plane of {Xsh, Qtt)- 


10 
























FIG. 12. CPC2: similar to Fig. 10 but in the plane of (g'*, g^^). 
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FIG. 13. CPVl: The ratio a/asM versus \^h for {g^^ g^) = (1, 0), (0,1), (l/-v/2, l/-v/2) 
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g®I gpl 



- 3 - 2-10 1 2 3 


FIG. 14. CPVl: Contours for the ratio (t/cjsm = 1 in the plane of (upper-left), 

{^ 3 H, Qt) (upper-right), and {g^,g^) (lower) for a few values of the third parameter. 
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LHC-14, Detector Level-ATLAS, CPV, ^ 13^=1 



LHC-14, Detector Level-ATLAS, CPV, Xg^sl 



FIG. 15. CPVl: Invariant mass distribution and for the decay products of the 

Higgs-boson pair. 


41 




















LHC-14, Detector Level-ATLAS, CPV, 



LHC-14, Detector Level-ATLAS, CPV, Xg^sl 



FIG. 16. CPVl: Angular distributions of and ARfji, between the two photons and between 

the two b quarks 
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aP, tP, 



FIG. 17. CPVl: The 25% and 50% sensitivity regions in the plane {gf, gf) bounded by three 
measurements of cross sections with basic cuts, > 2, and < 2 (the upper-left panel); 

with basic cuts, ARbh > 2, and ARi,b < 2 (the upper-right panel); and basic cuts, AR^^, AR^b > 2, 
and AR^^, ARbb < 2 (the lower panel). We assume that the measurements agree with the SM 
values with uncertainties of 25% and 50%, respectively. 
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25% & 50% forcr/o-sM = 1 
gst = l 

Detector Level -- ATLAS, LHC-14 

Basic Cuts - 

ARyy.bb > 2 

ARyy^ijb ^ 2 





-6 -4 -2 0 2 4 6 


FIG. 18. CPVl: similar to Fig. 17 but in the plane of {Xsh, 9t) 
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FIG. 19. CPVl: similar to Fig. 17 but in the plane of {X 3 H, 9t)- The left panels are for 


while those on the right are for gf = 0.5. 
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‘S 25% &50% fortr/tT-sM = 1 
Asb = l,gs^ = l.gP^ =0 
Detector Level - ATLAS, LHC-14 
Biisic Cuts - 

ARyy > 2 -- 

< 2 - 



25% &50% fortr/trsM = 1 

a = l,gst = l,gPt =0 

Detector Level - ATLAS, LHC-14 
Basic Cuts 

ARbb > 2 .— 

ARbb < 2 - 




FIG. 20. CPV2: The 25% and 50% sensitivity regions in the plane {g^^, g^) (with fixed X^h = 1, 
gf = 1, and gl^ = 0) bounded by three measurements of cross sections with basic cuts, ARy^ > 2, 
and AR^^ < 2 (the upper-left panel); with basic cuts, > 2, and AR^i^ < 2 (the upper-right 

panel); and basic cuts, AR^^, > 2, and AR^^, Ai?f,b < 2 (the lower panel). We assume that 

the measurements agree with the SM values with uncertainties of 25% and 50%, respectively. 
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FIG. 21. CPV2: similar to Fig. 20 but in the plane of {X^h, Qu) for = 1, gj; = g^^ = 0. 
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FIG. 22. 
9u = 0.5. 


CPV2: similar to Fig. 20 but in the plane of {Xsh, 9u) 9t 


1, 9t’ = 0, and 
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LHC-100, Detector Level-ATLAS 



LHC-100, Detector Level-ATLAS 



FIG. 23. Distributions in the invariant mass distribution and the angular separation 

in the decay products of the Higgs boson pair with detector simulation at the 100 TeV pp machine. 
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FIG. 24. CPCl at 100 TeV: Exactly the same as Fig. except that it is the case of the 100 
TeV pp machine. 
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